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In this article we focus on the parameterized complexity of the Multidimensional Binary Vector Assignment problem 
(called bMVA). An input of this problem is defined by m disjoint sets ,..., V^, each composed of n binary 

vectors of size p. An output is a set of n disjoint m-tuples of vectors, where each m-tuple is obtained by picking 
one vector from each set V’’. To each m-tuple we associate a p dimensional vector by applying the bit-wise AND 
operation on the m vectors of the tuple. The objective is to minimize the total number of zeros in these n vectors. 
bMVA can be seen as a variant of multidimensional matching where hyperedges are implicitly locally encoded via 
labels attached to vertices, but was originally introduced in the context of integrated circuit manufacturing. 

We provide for this problem FPT algorithms and negative results (ST//-based results, W[2]-hardness and a kernel 
lower bound) according to several parameters: the standard parameter k (i.e. the total number of zeros), as well as 
two parameters above some guaranteed values. 

Keywords: parameterized complexity, kernel, above guarantee parameterization, AND-cross composition, multidi¬ 
mensional binary vector assignment, wafer-to-wafer integration, locally encoded multidimensional matching 


1 Introduction 

1.1 Definition of the problem 

In this paper, we consider the parameterized version of the MULTIDIMENSIONAL Binary VECTOR AS¬ 
SIGNMENT problem (bMVA). An input of this problem is described by m sets ,..., V™, each of 

these sets containing n p-dimensional binary vectors. We note V® = {vl ,..., uG for all i G and 

for all j € [n] and r £ [p], we denote by u® [r] £ {0,1} the component of u]. 

In order to define the output of the problem, we need to introduce the notion of stack. A stack s = 

{vf, V 2 , ■ ■ ■, v^) is an m-tuple of vectors such that V* £ [m], uf £ V®. The output of bMVA is a set S of 
n stacks such that for all i,jG [m] x [n], Vj belongs to only one stack (in that case, the stacks are said 
disjoint). An example of an instance together with a solution is depicted in Figure || 

[m] stands for m}. 
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1/1 V2 1/3 S 


Figure 1: Example of bMVA instance with m = 3,n = 4,p = 6 and of a feasible solution S of cost c{S) = 17. 

We are now ready to define the objective function. We define the operator A that, given two p- 
dimensional vectors u and v, computes the vector w = A u[2] A u[2],..., u[p] A ^[p]). We 
associate to any stack s a unique vector Vg = Aig[m] ■ 

We define the cost of a binary vector v as the number of zeros in it. More formally, if v is p-dimensional, 
c{v) = p — X]r-e[p] We extend this definition to a set of stacks S = {si,... ,s„} as follows: 
c(‘5') = Sjgfn] Finally, the objective of bMVA is to obtain a set S of n disjoint stacks while 

minimizing c(S'). In the decision version of the problem, we are given an integer k, and we ask whether 
there exists a solution S of cost at most k. The problem is thus defined formally as follows: 
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Problem 1 Multidimensional Binary Vector Assignment (bMVA) 
Input: m sets of n binary p-dimensional vectors, an integer k 

Question: Is there a set S' of n disjoint stacks such that c(S) < k 1 


In order to avoid heavy notations throughout the paper, we will denote an instance of bMVA only by 
I[m, n,p, k], the notations of the sets and vectors being implicitly given as previously. 

1.2 Application and related work 

bMVA can be seen as a variant of multidimensional matching where hyperedges are implicitly locally 
encoded via labels attached to vertices. However, this kind of problem was originally introduced by Reda 


et al. (2009) in the context of semiconductor industry as the “yield maximization problem in wafer-to- 
wafer 3-D integration technology”. In this context, each vector Vj represents a wafer, which is seen 
as a string of bad dies (0) and good dies (1). Integrating two wafers corresponds to superimposing the 
two corresponding strings. In this operation, a position in the merged string is “good” when the two 
corresponding dies are good, and is “bad” otherwise. The objective of Wafer-to-Wafer Integration is to 
form n stacks, while maximizing their overall quality, or equivalently, minimizing the number of errors 
(depending on the objective function). In the following, we will denote by max^ 1-bMVA the dual 
version of bMVA where given the same input and ouput, the objective is to maximize np — c{S), the total 
number of ones. 

The results obtained so far concerning these problems mainly concern their approximability. It is 
proved in Dokka et al. (2013) that when to = 3, bMVA is NP-hard but admits a |-approximation. 
We can also mention Dokka et al. (2014) which provides a /(mj-approximation for general to, and an 
APX-hardness for to = 3. The main related article is Duvillie et al. (2015) where it is proved that 
maxJ2 1-bMVA has no 0{p^~'^) nor 0(TO^“'^)-approximation for any e > 0 unless P = NP (even 
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when n = 2), but admits a ^-approximation algorithm for any constant c G N, and is FPT when param¬ 
eterized by p (which also holds for bMVA). Notice that one of the reductions provided by puvillie et al. 
( 2015| ) is a parameter-preserving reduction from the CLIQUE problem to max 1-bMVA, immediately 
proving W[l]-hardness for max Y 1-bMVA when parameterized by the objective function. This is why 
our motivation in this paper is to consider the parameterized complexity of bMVA. As we will see in the 
next section, we provide an analysis for several parameters related to this problem. 


For fo rmal definitions and detailed concepts on Fixed-Parameter Tractability, we refer to the mono¬ 
graph of Downey and Fellows (2013). Moreover, in order to define lower bounds on the running time of 


parameterized algorithms, we will rely on the Exponential Time Hypothesis (ETH) of [Impagliazzo et al 


( 2001 ), stating that 3 -Sat cannot be solved in where n is the number of variables (0*(.) hides 

polynomi al terms). For more result s about lower bounds obtained under ETH, we refer the reader to the 
survey of pokshtanov et ah (2011a). 


1.3 Parameterizations 


One of the main purposes of Fixed-Parameter Tractability is to obtain efficient algorithms when the con¬ 
sidered parameter is small in practice. When dealing with the decision version of an optimization problem, 
the most natural parameter is perhaps the value of the desired solution (e.g. k for bM VA). Such a param¬ 
eter is often referred to in the literature as the “standard parameter” of the problem. In some cases, this 
parameter might not be very interesting, either because it usually takes high values in practice, or because 
FPT algorithms with respect to this parameterization are trivial to find. When this happens, it is possible 
to obtain more interesting results by subtracting to the objective function a known lower bound of it. For 
instance, if one can prove that any solution of a given minimization problem is of cost at least B, then one 
can ask for a solution of cost B + c and parameterize by c. This idea, called “above guarantee parame¬ 
terization” was introduced by Mahajan and Raman (1999) and first applied to Max Sat and Max Cut 
problems. It then became a fruitful line of research with similar results obtained for many other problems 
(among others, see jCygan et al. (2012); Gutin et al. (2006); Gutin and Yeo (2012); Mahajan et al. (2009)). 


In this paper, we analyze the parameterized complexity of bMVA using three types of parameters. The 
first one is the standard parameter k: the number of zeros to minimize in the optimization version of the 
problem. Then, three natural structural parameters: m, the number of sets of the input, n, the number of 
vectors in each set, and p, the size of each vector. The last two parameters are above guarantee parameters. 


As we said previously, we alreadyproved in Duvillie et al. (2015) that bMVA is FPT parameterized 


by p. As we will notice in Lemma g that we can obtain p < k after a polynomial pre-processing step, 
this implies that bMVA is also FPT with its standard parameter. Our idea here is to use this previous 
inequality in order to obtain smaller parameters. Thus, we define our first above guarantee parameter 

(p = k-p. 


Finally, in order to define our last parameter, we first need to describe the corresponding lower bound 
B, that will represent the maximum, over all sets of vectors, of the total number of zeros for each set. 
More formally, we define B = maxjgj^j c(F®) where c(F*) = Y'j=i Since we perform a bit-wise 
AND over each m-tuple, it is easily seen that any solution will be of cost at least B. Thus, we define our 
last parameter Qij = k — B. 
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1.4 Our results 

In the next section, we present some pre-processing rules leading to a kernel of size 0(k^m), and prove 
that even when m = 3 we cannot improve it to unless NP C coNP/poly (remember that bMVA 
is known to be FPT when parameterized by p). Section 0 is mainly focused on results associated with 
parameter C,s'- we prove that bMVA can be solved in log(n)^^ while it is W[2]-hard when parame¬ 
terized by C,B only, and cannot be solved in iog(")) nor in assuming ETH. In 

Section Q, we focus on the parameterization by C,p: we show that when n = 2, the problem can be solved 
in single exponential time with this parameter, but is not in XP for any fixed n > 3 (unless P = NP). 
The reduction we use also shows that for fixed n S N, the problem cannot be solved in 2°^^^ (and thus 
in 2°^“’®^) unless ETH fails, which matches the upper bound obtained in Section || A summary of our 
results is depicted in the following table. 


Positive results 

Negative results 

0{k^m) kernel (Thm. 1) 

no kernel for m = 3 unless NP C coNP/poly (Thm. 3) 

0 *( 4 Ceiog(n)) algorithm 
(Thm.|) 

W [2]-hard for only (Thm. ^ 

no 2°(?e) nor under ETH (Thm. |) 

no 2°(^) for fixed n under ETH (Thm. | 8 |) 

0*{d^p) algorithm 
for n = 2 (Thm. Q) 

NP-hard for Cp = 0 and fixed n > 3 (Thm.|^ 


This article is the complete version of Bougeret et al. (2015) where some proofs were omitted due to 
space limitations. 


2 First remarks and kernels 


Let us start with two simple lemmas allowing us to bound the size of the input. Notice first that it is 
not always safe to create a 1-stack (i.e. a stack with ones on every component) when possible. Indeed, 
in instance 1/1 = {(111), (101), (Oil)}, = {(111), (101), (110)}, = {(111), (Oil), (110)}, 

depicted in Figure ^ no optimal solution creates a 1-stack. However, as we will see in Lemma[|, creating 
1 -stacks becomes safe if n > k. 




L2 


y3 


Opt 



c(VsJ = 1 
c(Vs^) = 1 
c(Vss) = 1 


Figure 2: Example of bMVA instance such that no optimal solution creates a 1-stack. 


Lemma 1. There exists a polynomial algorithm which, given any instance 

I[m, n,p, k] of bMVA, either detects that X is a negative instance, or outputs an equivalent instance 
I'[m,n',p,k] such thatn' < k. 
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Proof: Let I[m, n,p,k] be an instance of bM VA, and suppose that n > k. Let us write a polynomial pre¬ 
processing rule that either detects that I is a no instance, or compute an equivalent instance X" [m,n",p, k] 
with n" = n — 1. 

Notice first that there exists at least a 1-vector in every set L®. If not, I is a no instance as any solution 
would be of cost at least n > k. \t is now safe to create a 1-stack, obtaining a remaining instance I" 
with n" = n — 1. Indeed, if I is a yes instance, then there must exist at least one 1-stack in the solution 
(otherwise the cost would be at least n > k), and thus the remaining instance is also a yes instance. As 
the converse is trivially true, the rule is safe. Applying it at most n — k times finally leads to the desired 
upper bound. □ 


Lemma 2. There exists a polynomial algorithm which, given any instance 

I[m,n,p, k] of bMVA, either detects that X is a negative instance, or outputs an equivalent instance 
X'[m, n,p', k] such that p' < k. 


Proof: Let I[m, n, p, k] be an instance of the problem, and suppose that there exists r € [p] such that for 
all {i,j) € [m] X [n] we have u][r] = 1. In other words, the component of all vectors of all sets is 
a 1. In this case, it is clear that all vectors of any set of n stacks obtained from X will also contain a 1 
at the component. Hence, we can modify X[m, n,p, k] into I'[to, n,p', k] with p' < phy dropping 
all such components for all vectors. It is clear that this rule is safe since the cost of any solution remains 
unchanged, and it can be applied in polynomial time. After applying this rule, for all r € [p'] there exists 
{i, j) € [to] X [n] such that u* [r] = 0. This immediately implies that the cost of any solution is at least p', 
and thus if p' > fc the algorithm detected a no instance. □ 


Given the two previous lemmas, we can suppose from now on that for any instance of bMVA we have 
n < k and p < k. This immediately implies a polynomial kernel parameterized by k and to. 

Theorem 1. bMVA admits a kernel with 0{kfm) bits. 

Let us now turn to the main result of this section. To complement Theorem we show that even 
when TO = 3, we cannot obtain a polynomial kernel with the smaller parameter p under some classical 
complexity assumptions (notice however that the existence of a polynomial kernel in k only is still open 
when TO is not fixed). Notice also that as bMVA was known to be FPT when parameterized by p, proved 
in Duvillie et al. (2015), it was a natural question to ask for a polynomial kernel. 

In order to establish kernel lower bounds, we use the concept of AND-cross-composition of Bodlaender 
et al. (2014), together with the recently proved AND-conjecture of prucka (2015). In the following, a 
parameterized problem is a subset of E* x N, where S* is the set of words over some finite alphabet S. 
Definition 1 (Polynomial equivalence relation according to Bodlaender et al. (2014)). An equivalence 
relation TZ on E* is called a polynomial equivalence relation if the two following conditions hold; 


• There is an algorithm that given two strings x,y G E*, decides whether x and y belong to the same 
equivalence class in ([a: I + time. 


• For any finite set S C E*, the equivalence relation TZ partitions the elements of S into at most 

{maXx^s\x\)^^^'’ classes. 


Definition 2 (AND-cross-composition according to Bodlaender et al. (2014)). Let L C E* be a set and 
let Q C E* X N be a parameterized problem. We say that L AND-cross-composes into Q if there 



















6 


M. Bougeret, G. Duvillie, R. Giroudeau, R. Watrigant 


is a polynomial equivalence relation TZ and an algorithm which, given t strings belonging to the same 

equivalence class of TZ, computes an instance (a;*, k*) G E* x N in time polynomial in 

that; 


• {x*, k*) G Q a;i G L for alH G {1 ,..., f} 

• k* is bounded by a polynomial in max\^i | + log t 


Theorem 2 ( pruckeii ( l2015| )). If some set L C E* is NP-hanf and L AND-cross-composes into a pa¬ 
rameterized problem Q, then there is no polynomial kernel for Q unless NP C coNP/poly. 


Theorem 3. Even for m = 3, bMVA parameterized by p does not admit a polynomial kernel unless 

NP C coNP/poly. 


Proof: The proof is an AND-cross-composition inspired by the NP-hardness reduction for bMVA pro¬ 
vided by Dokka et al. (2013). More precisely, we cross-compose from a sequence of instances of 3- 
DIMENSIONAL PERFECT MATCHING. According to Dehnition ^ and Theorem ^ this will imply the 
desired result. 3-DlMENSlONAL PERFECT MATCHING is formally defined as follows: 


Problem 2 3-Dimensional Perfect Matching 


Input; Three sets X, Y and Z of size n, a set of hyperedges S C X x Y x Z 

Question: Does there exist a subset S' f S such that: 


• for all e, e' G S' with e = {x,y,z) and e' = {x',y',z'), we have x f x', 
y f y' and z z' (that is. S' is a matching) 

• |S"| = n (that is. S' is perfect) 


Let {Xi, Yi, Zi, Si),..., {Xt,Yt, Zt, St) be a sequence of t equivalent instances of 3-DPM, with re¬ 
spect to the following polynomial equivalence relation: {X, Y, Z, S) and {X', Y',Z', S') are equivalent 
if I AT I = \X'\ (and thus |E| = \Y'\ = \Z\ = \Z'\ = |Ar|), and \S\ = |S"|. In the following we denote 
by n the cardinality of the sets Xi (and equivalently the sets Yi and Zi), and by m the cardinality of the 
sets S^. Moreover, for alH G t} we dehne Xi = {xip,Y^ = {yi^i,yi^n}, Zi = 

{zi^i,..., Zi^n}, and Si = ..., Si^m}- We also assume that t = 2^ for some g G N (if it is not the 

case, we add a sufficiently number of dummy yes-instances). 

In the following we construct three sets {X*, Y*,Z*) of nt vectors each: X* = {x* Y* = 

{y*jViZi’ T where each vector is composed of p* = m -f 2mq components. 

Let us hrst describe the first m components of each vector. For alH G {1,... ,t}, j G {1,..., n} and 
fc G {1,..., m} we set: 




[k] 


1 if the hyperedge Si^k contains Xij 
0 otherwise 
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ViAk] 


1 if the hyperedge Si^k contains yij 
0 otherwise 


,[k] 


1 if the hyperedge Si^k contains Zij 
0 otherwise 


Then, for alH G {1,..., f}, we append two vectors bi and bi to all vectors {x* {y* 

and {z* The vector bi is composed of mq coordinates, and is defined as the binary representa¬ 

tion of the integer i, where each bit is duplicated m times. Finally, bi is obtained by taking the complement 
of bi (i.e. replacing all zeros by ones, and conversely) as depicted in Figure |^. It is now clear that each 
vector X* ^ (resp. z* j) is composed of p* = to + 2mq coordinates. Thus, the parameter of the input 
instance is a polynomial in n, m and log t whereas the total size of the instance is a polynomial in the size 
of the sequence of inputs, as required in cross-compositions. It now remains to prove that {X*,Y*, Z*) 
contains an assignment of cost k* = nt{mq -f to — 1) if and only if for all i G {1 ,... ,t}, Si contains a 
perfect matching S''. 

• Suppose that for all i G {1, ■ ■ ■ ,t} we have a perfect matching S' C Si. W.l.o.g. suppose that 

S' = ..., Then, for each j G {l,...,n}, we have = {xij^,yij 2 , Zijs) for some 

ji) j 2 , ja G n}. We assign with y*j^ and z*^. It is easy to see that the cost of this 

triple is TO — 1 -f mq. Indeed, they all have a one at the coordinate, corresponding to the 
hyperedge of Si (and this is the only shared one, since we can suppose that all hyperedges are 
pairwise distinct), and they all contain the same vectors bi and bi. Summing up for all instances, we 
get the desired solution value. 

• Conversely, first remark that in any assignment, the cost of every triple , z*^ ) is 

at least to — 1 -f mq, and let us prove that this bound is tight when (1) all elements are chosen 
within the same instance, i.e. ii = i 2 = is = F and (2) this triple corresponds to an element 
of Si, i.e. {xi,ji,yi,j 2 : Zij^) G Si. Indeed, suppose first that ii ^ 12 . Then, since the binary 
representation of zi and 12 differs on at least one bit, it is clear that the resulting vector is of cost at 
least m{q -f 1) > mq + m — 1. Now if ii = 12 = = F then the result is straightforward, since at 

most one hyperedge of Si can contain y*^ and z*^ . Finally, using the same arguments as 

previously, we can easily deduce a perfect matching S’' C Si for each i G {1,..., f}, and the result 
follows. 


□ 


3 Parameterizing according to (b 

In this section, we present an FPT algorithm when parameterized by and n (recall that both Cs and 
n are smaller parameters than the standard one k, since k = B + Cb and n < k in any reduced instance). 
Notice first that it is easy to get an C)*( 2 ‘^e(iog(n)-i-iog(p))^ algorithm. Indeed, by considering a set i G [to] 
where c(C') = B, and guessing the positions of the (s new zeros (among np possible positions) that 
will appear in an optimal solution, we can actually guess in 0*{{np)^’^) the vectors {us*} of an optimal 


M. Bougeret, G. Duvillie, R. Giroudeau, R. Watrigant 
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Figure 3: Example of Cross Composition construction from four equivalent instances of 3 -Dimensional 
Perfect Matching {Xi,Yi, Zi, yi,i, 21,3), t/1,2, 21,2), (*1,2, yi,3, 21.3), (*1,3, t/1,3, 21,3)}), 

(X2,12, Z2, {(a;2,i, y2,i, 22,1), (*2,2, j/2,2,22,2), (*2,3, t/2,3, 22,3), (2:2,3, 2/2,2, 22,1)}), 

(^3, Ys, Z3, {(a; 3 ,i, t/3,1,23,1), (*3,2, j/3,1,23,1), (*3,3, 2 / 3 ,2, 23,2), (2:3,3, 2 / 3 , 3 , 23,3)}) 

(X4, > 4 , ^4, {(2:4,1, 2 / 4 , 2 , 24,2), (2:4,2,2/4,1,24.1), (2:4.2,2/4,2, 24,2), (2:4,3, 2 / 4 . 3 , 24.3)}) into an instance of bMVA 

with m = 3 , n = 12 and p = 20 . 
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solution, and it remains to check in polynomial time that every can be “matched” to {ug*}. Now we 
show how to get rid of the log(p) term in the exponent. 

Theorem 4. bMVA can be solved in ^og{n)y 

Proof: Let X[m,n,p, k] be an instance of our problem and, w.l.o.g., suppose that is a set whose 
number of zeros reaches the upper bound B, i.e. c{V^) = B. The algorithm consists in constructing a 
solution by finding an optimal assignment between and ,..., 14™, successively. 

We first claim that we can decide in polynomial time whether there is an assignment between and 
which does not create any additional zero. 

To that end, we create a bipartite graph G with bipartization (A, B), A = {oi,..., a„}, B — {bi,..., 
bn}, and link and bj^ for all (ji, J 2 ) G [n-] x [n] iff assigning vector from and vector Vj^ from 
does not create any additional zero in (vj^ A = vj^). 

If a perfect maching can be found in G, then we can safely delete the set 14^ and continue. In order to 
avoid heavy notations, we consider this first step as a polynomial pre-processing, and we re-label 14* into 
14*“^ for all i G {3,..., m} (and m is implicitly decreased by one). 

In the following, we suppose that the previous pre-processing step cannot apply (i.e. there is no perfect 
matching in G). Intuitively, in this case any assignment (including an optimal one) between and 
must lead to at least one additional zero in V^. In this case, we perform a branching to guess one couple 
of vectors from x 14^ which will induce such an additional zero. More formally, we branch on every 
couple (ji, J 2 ) G [n] x [n], and create a new instance as a copy of I in which is replaced by A ■ 
This operation increases c(V^) by at least one, and thus B by at least one as well. If we denote by I' this 
new instance, we can check that a solution of cost at most k for I' will immediately imply a solution of 
cost at most k for I, as I' is constructed from I by adding some zeros. The converse is also true as one 
assignment we enumerate corresponds to one from an optimal solution. 

As the value of B in this branching increases by at least one while we still look for a solution of cost k, 
this implies that this branching will be applied at most Cs times. Summing up, we have one polynomial 
pre-processing and one branching of size n? which will be applied at most (b times. The total running 
time of this algorithm is thus bounded by O* (4^® iog(")), 

□ 

Despite its simplicity, we now show that, when considering each parameter (n and ((g) separately, this 
algorithm is the best we can hope for (whereas the existence of an O* (2^) algorithm is still open). Indeed, 
we first show in Theorem^ that the linear dependence in and log{n) in the exponent is necessary 
(unless ETH fails), and also that we cannot hope for an FPT algorithm parameterized by only unless 
FPT = W[2] (Theorem |]). Finally, as we will see in the next section (Theorem ||), this result is matched 
by a lower bound when n G N is fixed. We now present a reduction from the HITTING Set 

problem which produces an instance of bMVA. 


Problem 3 Hitting Set 


Input: ^ subsets i?i,..., Rm of [n], and an integer k 

Question: Is there a set Rofk elements of [n] such that RC\ Ri for any i G [m] 1 
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Figure 4: Example of reduction from an instance of Hitting Set consisting of four subsets of [n = 7 ] : i?i = 
{ 1 , 2, 7 } , R2 = { 1 , 3 , 4, 5 } , -R3 = {2, 4, 5 } , -R4 = { 5 , 6, 7 }, and an integer fc = 2, to an instance of bMVA with 
m! = 7>,n' = n,p — n. 


Lemma 3. There is a polynomial reduction from HITTING Set to bMVA that given an instance com¬ 
posed of m subsets of [n] and an integer k, constructs an instance of bMVA X[m', n',p', k'] such that 
n' = n and C,b = k. 


Proof: Let Ri,Rm be subsets of [n], and k G N. We construct m sets 17^,of n vectors each, 
where, for all i G [m] we have 17* = {v\, each vector being composed of n components. For all 

i G [m] and all j G [n], if j G Ri, then the vector u* is composed of ones everywhere except at the 
component. If j ^ Ri, then u* is a 0-vector (i.e. a vector with zero in every component). We also add a 
set 17* composed of {n — 1) 0-vectors and one 1-vector as depicted in Figure^. 

For this constructed instance, it is clear that B = n{n — 1) because of the set 17*. In other words, any 
assignment will lead to a solution with (n — 1) 0-vectors, and thus with at least n{n — 1) zeros. We will 
actually show that this instance has a solution with n(n — 1) -f fc zeros if and only if i?i,..., Rm has a 
hitting set of size k. By the foregoing, we only need to focus on the only vector of each set which is 
assigned to the 1-vector of 17*. 

Let J C [n] be a hitting set of size k. By the definition of a hitting set, for all i G \m\, there exists 
ji G J r\ Ri- Thus, for all i G [m], we select the vector u*. from the set 17* to be assigned to the 1- 
vector of 17*. By construction, this vector has only one zero at the jf^ component, which implies that the 
conjunction of all such vectors /\i^i will have a 1 everywhere except at the components corresponding 
to J. We thus have the desired number of zeros in our solution. 

4= Conversely, for each i G [m], let ji G [n] be the vector from 17* which is assigned to the 1-vector of 
17*. Since the resulting conjunction of all these vectors has only k zeros, v^. cannot be a 0-vector, and 
we thus have ji G Ri- Using the same arguments as previously, {uji}ig[m] corresponds to a hitting set of 
i?i,..., Rm of size k. 

As we seen previously, B = n{n — 1) for the obtained instance, k' = n(n — 1) -\- k, (which implies 
Cb = k), and the size of all sets is n, as desired. □ 


As we ca n see, the reduction is parame ter-preserving for From the W [2]-hardness of HITTING Set 
provided by Downey and Fello^ (2013), we have the following; 


Theorem 5. bMVA is 'W^]-hard when parameterized by (jg. 
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As said previously, we also use this reduction to show the following result: 

Theorem 6. bMWA cannot be solved in nor unless ETH fails. 

Proof: To show this, we u se a the previous hardness r esult, but using a constrained version of the HITTING 
Set problem obtained by Lokshtanov et al. ( 201 lb ), where the element set is [fc] x [fc] and can thus be 
seen as a table with k rows and k columns: 


Problem Ak xk Hitting Set 


Input; 

Question: 


An integer k, and Ri, ...,Rt C [A:] x [k] 

Is there a set R containing exactly one element from each row such that RPi Ri f 
for any i G [f] ? 


Lokshtanov et al. ( 2011b ) show that assuming ETH this problem cannot be solved in 
(whereas a simple brute force solves it in O* Notice that we can modify the question of this 

problem by dropping the constraint that S contains at least one element from each row. Indeed, let us add 
to the instance a set of k sets ..., i?^}, where i?' contains all elements of row i for i G [fc]. Now, 
finding a (classical) hitting set of size k on this modified instance is equivalent to finding a solution of size 
k for the original instance of fc x fc-HiTTiNG Set. Moreover, it is easy to check that a 
algorithm for this relaxed problem would also contradict ETH. To summarize, we know that unless 
ETH fails, there is no 2°^^ log(fe))^.^ 0 ( 1 ) algorithm for the classical HITTING Set problem, even when the 
ground set has size This allows us to perform the reduction of Lemma on these special instances, 
leading to an instance I[m', n',p', fc'] with associated parameter Cg such that Cb = k and n' = k^. 
Suppose now that there exists an algorithm for bMVA running in 2°^^®^ iog(") 77 J _|_ _|_p) 0 (i) pj^ing 
the reduction above, we would be able to solve the instance of A: x A;-HlTTlNG Set in 2°^^^ iog(fc^)^o(i)^ 
and thus in which would violate ETH. A similar idea also rules out any algorithm 

running in under ETH. □ 


4 Parameterizing accorcding to Cp 

We now consider the problem parameterized by = A; — p (recall that p < k). Notice that one motivation 
of this parameterization is the previous reduction of LemmaH from HITTING Set. Indeed, when applied 
for n = 2, it reduces an instance of VERTEX COVER to an instance of bMVA with k = p + Cp where 
(p is equal to the size of the vertex cover. Our intuition is confirmed by the following result: we show 
that when parameterized by (p, the problem is indeed FPT when n = 2 (Theorem L). We complement 
this by showing that for any n > 3, it becomes NP-hard when (p = 0 (Theorem K), and is thus even 
not in XP. The reduction we use even proves that for any fixed n > 3, the problem cannot be solved in 
2 o(fc) (and thus in 2°^^®)) unless ETH fails, while the algorithm of Theorem^ runs in In 

the following, n-BMVA denotes the problem bMVA where the size of all sets is fixed to some constant 
n G N. 

4 .1 Positive result torn = 2 

In this subsection, we prove that 2 -bMVA is FPT parameterized by ^p. To do so, we reduce to the Odd 
Cycle Transversal problem (OCT for short). In this problem, given a graph G = {V, E) and an 
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integer c G N, the objective is to decide whether there exists a partition {X, 81 , 82 ) of V with jXl < c 
such that 81 and 82 are independent sets. 

We hrst introduce a generalized version of OCT, called BIP-OCT. In this problem, we are given a set of 
vertices V, an integer c, and a set of m pairs {Ai, Bi),(Am, Bm) with Ai,Bi C V for all i G [m] and 
AiCi Bi = 0. Informally, each pair (Ai, Bi) can be seen as a complete bipartite subgraph. The output of 
BIP-OCT is described by a partition (X, 81 , 82 ) of V such that for any i G [m], either (Ai \ X G 81 and 
Bi\X G 82 ) or (Ai \ X G 82 and Bi\X G 81 ). The question is whether there exists such a partition 
with | 2 f| < c. As we can see, if all Ai and Bi are singletons (and thus form edges), then BIP-OCT 
corresponds to OCT. Notice that in the following, the considered parameter of OCT and BIP-OCT will 
always be the standard parameter, i.e. c. We hrst show that there is a linear parameter-preserving reduction 
from 2 -bMVA parameterized by Cp to BIP-OCT, and then that there is also a linear parameter-preserving 
ti'ansformation from BIP-OCT to OCT. 

Lemma 4. There is a linear parameter-preserving reduction from 2-bMVA parameterized by Cp to BIP- 
OCT. 

Proof: Let I[m,2,p,p -f (p] be an instance of 2 -bMVA (i.e. in which every set contains only two 
vectors), and let us construct an instance X' of BIP-OCT, such that X has a solution of cost p -f Cp iff 
has a solution of size Cp. Notice hrst that we can suppose that for any i G [m] and any r G [p], we cannot 
have both vl [r] = 0 and V 2 [r] = 0 as otherwise any stack s from any solution would have Vg [r] = 0 , and 
thus we could safely remove such a component r from the instance (and decrease k and p by one). 

Let the vertex set of X' be [p]. Then, for all i G [m], let us dehne Ai = {r|uj[r] = 0}, and Bi = 
{r|u 2 [r] = 0} as depicted in Figure By the foregoing, and as required in an instance of bip-OCT, we 
have AiG Bi = 0. Let us prove that X has a solution of cost p -f Cp iff B has a solution (X, 81 , 82 ) with 

1^1 < Cp- 

Let 8 = {si,S 2 } be a solution of X of cost p -f C.p. Let X = jrlrisj?"] = Us 2 [r] = 0}, 81 = 

{r|r;si H = 0 Vg.^ [r] = 1}, and 82 = {r|usi W = 1 Vg.^ [r] = 0}. Notice that (X, 81 , 82 ) forms 

a partition of [p] (as we cannot have a tq with Vg^^ [pq] = Vg.^ [pq] = 1 , as this would imply that all the nm 
vectors have p[po] = 1, and such cooordinates have been removed from the instance in Lemma and 
\X\ = Cp. It remains to prove that (2f, S' 2 ) is a feasible solution of X'. Let i G [m]. Without loss of 

generality, let us suppose that v\ has been added to si and v\ b^s been added to S 2 . Let r G Ai\X. Since 
r G Ai, we have v\ [r] = 0, and thus [p] = 0. Since r ^ X, we have Vg^ [p] = 1. Thus, r G 81 , which 
proves Ai \ C 81 . Similarly, we can prove that Bi \ X G 82 - 

Let (X, 81 , 82 ) be a solution of X' with |2f| < Cp- Let si be such that Vg^ [p] = 0 iff p G X or 

r G 81 , and let S 2 be such that Vg.^ [p] = 0 iff p G AT or p G 52 . It remains to prove that the solution 

8 = {si, S 2 } is feasible, which immediately implies that its cost is p -f Cp- Let i G [m]. Without loss of 
generality, let us suppose that \ X C 5i and Bi\X G 82 . We now claim that v\ can be assigned to 
Si and V 2 can be assigned to S 2 without creating any new zero. To do so, let us show that for all p G [p], 
we have v\[r] = 0 => PsJp] = 0 (resp. v^r] = 0 => = b). Indeed, let p G [p] such that 

v\[r] = 0. Then by construction, it means that p G Ai. Thus, by dehnition of the solution (X, 81 , 82 ), 
it means that either p G X or p G 5i, which implies Vg^ [p] = 0 as desired. Similar arguments show that 
p® [p] = 0 => Ps 2 [p] = 0 for all p G [p]. □ 


Lemma 5. There is a linear parameterized reduction from BIP-OCT to OCT. 
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001111011 


111010111 


110100110 



011101111 


110101110 


111111000 


001001000 


Figure 5: Example of reduction from an instance of 2 -bMVA, with n = 2 , m = 3 ,p = 9 , admitting a solution 
of cost p + 2 , to an instance of BIP-OCT, with \V\ = p, Ai = { 3 , 6} , Bi = { 1 , 5 } , A2 = { 2 , 3 , 7 } , -B2 = 
{ 3 , 5 , 9 } , A3 = { 4 , 6} , -63 = { 7 , 8, 9 }, admitting the partition ({ 5 , 9 } , { 3 , 4 , 6} , { 1 , 2 , 7 , 8}) as solution of cost 2 . 
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Proof: Let I = be an instance of BIP-OCT. Let us construct a graph G" = {V',E') 

which contains an odd cycle transversal of size c if and only if I has a solution of size c for BIP-OCT. 
Observe first that we cannot simply set V' = V and E' = lJi6[m] asA- Indeed, if for example 

Ai = {2,..., n}, Bi = {1}, A 2 = {2, ..., j} and B 2 = n}, defining G' as above would 

lead to an odd cycle transversal of size one, as removing only vertex { 1 } makes the graph bipartite with 
bipartization (A 2 , i? 2 )- However, this solution is not feasible for BIP-OCT as \X = Ai, and Ai % A 2 
and Ai % B 2 . Intuitively, we have to prevent solutions of G' from splitting sets Ai\X (and Bi \ X) 
between the two parts of the bipartization. To do so, we will construct G' as described above, and then 
we "augment" each bipartite graph by adding c+1 new vertices on each side. More formally, we start by 
setting y' = L as said before, and for all i € [m], we create two sets of c -f 1 new vertices A', B[. We 
then set E' = U*eH,aGA,uA'.faeB,uB'{a> 

Let us now prove that X contains a solution of size c for BIP-OCT if and only if G' contains an odd 
cycle transversal of size c. 

=> Let {X, 81 ^ 82 ) be an optimal solution of I. We define a partial solution X', 8[,82 of G' by setting 
X' = X and 8 [ = 81 for I S {1,2} (the solution is partial in the sense that it remains to assign vertices 
of A' U S', for all i € [m]). Let i € [m]. If \ X = 0 and Bi\ X = 0, then we add (arbitrarily) 
A' to 8 'i and to 82 . Otherwise, if \ Jf 7 ^ 0 and is added to 81 , we add A' to S'} and to S}/ 
with I, I' G { 1 , 2 }, V 7 ^ /, and if \ X 7 ^ 0 and is added to 81 , we add to S} and A' to S}/ with 
1 , 1 'G { 1 , 2 }, I'^l. 

This new solution has the same size (|X'| = |X|) and we claim that it is an odd cycle transversal of G'. 
Indeed, let us check that any edge [u, w} G E' such that {u, v} fl X' = 0 is not entirely contained in a 
S}. If {tt, u} is an edge of a complete bipartite of X, i.e. if there exists i such that u G Ai and v G Bi, 
then by definition of the solution {X, 81 , 82 ) it is straightforward that u and v are not both in S} nor in 
82 - Otherwise, if {tt, u} is adjacent to one or two of the new vertices, let i be such that u G A'. If u G S', 
then the solution is valid as A' and B' are never added to the same set S'}, I G {1, 2}. Otherwise, we 
necessarily have v G Bi. Let I G {1, 2} be such that Bi\ X (which is not empty) has been added to S}. 
In this case Ai (and thus u) has been added to S},, with V ^ 1. 

Let {X', 8 'i, 82 ) be an optimal solution of G'. For any i G [m], let Ai = {Ai U A') \ X' and 
Bi = {Bi U S') \ X'. A first observation is that Ai ^ % and 7 ^ 0 as \Ai U A'| = \Bi U i3'| > c and 
I AT'I < c. A second observation is that for any u and v G Ai, u and v are in the same set S} for some 
I G {1, 2}. Indeed, suppose by contradiction that u G 8 'i and v G 82 - As Bi 7 ^ 0, there exists b G Bi and 
I G {1, 2} such that b G 8 'i. As all the edges of the complete bipartite subgraph on {Ai U A', Bi U B'i) 
belong to E' , we have {u, h} G E' and {u, b} G E' , and thus S'} contains both endpoints of an edge of 
E', which is a contradiction. In the same way, we can prove that for any i G [m], and any u and v G Bi, 
u and V are in the same set S} for some I G {1,2}. 

Thus, according to the two previous observations, for any i G \m\ we can define A^. G {1,2} and 
\g. G {1,2} such that Ai C and Bi C , with A^. ^ 

Let us now define X = X' n L, Si = S} fl L, and S 2 = S 2 n V, and check that this is a valid solution 
of X. Let i G [m]. Observe first that \ X C Ai, and thus either A^ \ X = 0, or A^ \ X C 8 \^ . As the 
same fact also holds for Bi \ X, and as A^. 7 ^ A^., the constraint {Ai \ X C Si and Bi \ X C 82 ) or 
(Ai \ X C 82 and Bi\X C Si) is respected, and the solution is feasible, which concludes the proof. 

□ 



Multidimensional Binary Vector Assignment problem 


15 


As Odd Cycle Transversal can be solved in O*(2.3146'^), proved by f^okshtanov et al. (2014), 
and since our parameters are exactly preserved in our two reductions, we obtain the following result: 


Theorem 7. 2 -bMVA can be solved in 0*{d‘^r) where d < 2.3146 is such that OCT can be solved in 
0 *{d^). 


4.2 Negative results forn>3 

We now complement the previous result by proving that the problem is intractable with respect to the 
parameter Cp for larger values of n. 

Theorem 8. For any fixed n > 3, n-BMVA is not in XP when parameterized by (unless P = NPj, 
and cannot be solved in (unless ETH fails). 


Proof: Let x ^ 3. We present a reduction from x-COLORiNG. Given a graph G = (F, E), this problem 
consists to ask for a mapping / : V —[x] such that for all {m, t;} G E we have f{u) f(.v). Let 
E = {ei, ..., 6 ^ 0 } V = [nc]. Let us construct an instance I of n-BMVA with n = x, P = no, 

m = me and such that G admits a x-coloring iff I has a solution of cost p (i.e. (p = 0). To each edge 

Ci = {m, v} G E,i G [ma], we associate a set F* with |F*| = x. where: 

• v\ represents the vertex u, that is v\ [n] = 0 and v\ [r] = 1 for any r G [no], r u, 

• V 2 represents the vertex v, that is v^fv] = 0 and v\\r\ = 1 for any r G [na], r v, 


for all j G {3,..., x}, f* is a 1-vector, i.e. it has a 1 at every component. 


An example of this construcion is depicted in Figure || Let us now prove that G admits a x-coloring iff I 
has a solution of cost p = nc. 

=> Let Sj G V, j G [x] be the x color classes (notice that the Sj are pairwise disjoint, some of them 
may be empty, and Uj 6 [x] ~ each Sj we associate a stack sj such that [r] = 0 iff r C Sj. 

It remains to prove that the solution S = {si,..., is feasible, as its cost is exactly p by construction. 
Let us consider a set F* where v\ (resp. vlf) represents a vertex u (resp. v). As {n, u} is an edge of G, 
we know that u and v have two different colors, i.e. that u G Sj and v G Sj/, for some j,j' G [x] with 
j j'. Thus, we can add v\ to stack Sj, v\ to stack Sj', and the x “ 2 other u* (j > 3) in an arbitrary 
way. Since the only 0 in v\ (resp. V 2 ) is at the (resp. component, we have v\ A Vg^ = Vg^ (resp. 
V 2 A Vg^, = Vg^,), which proves that S is feasible. 

<t= Let S = {si,..., S;^} be the stacks of an optimal solution. For j G [x], let Sj = {r G [pJIt's [f] = 

Ji 

of F. Moreover, as for any i G [m], v\ and v. 
ing vertices have been assigned to different colors, and thus each Sj induces an independent set, which 
completes the reduction. 

It is known, thanks to Impagliazzo et al. ( 2001 ) that there is no 2°ll^b algorithm for deciding whether 
a graph G = (F, E) admits a x-COLORiNG, for any x ^ 3 (under ETEl). As we can see, the value of 
the optimal solution for n-BMVA in the previous reduction equals the number of vertices in the instance 
of x-COLORiNG, which proves that n-BMVA cannot be solved in 2°1^1 for any n > 3. 

□ 


0}. Notice that lj?=i ‘b'i = V, and as I is of cost p, all the Sj are pairwise disjoints and form a partition 

\ and V 2 have been assigned to different stacks, the correspond- 
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yl.2 yl,3 yl,4 y2,3 y3,4 y2,5 y4,5 g 



Figure 6: Example of reduction from a positive instance of ^-COLORING, with x = 3 , V = [ 5 ],-B = 
{{1, 2} , {1, 3 } , {1, 4} , {2, 3 } , { 3 , 4} , {2, 5 } , {4, 5 }}, to an instance of bMVA with m = |i5| = 6, n = x = 
3 , p = |1^1 = 5 admitting a solution of cost p. 

Finally, remark that as for the parameterization by p, one could ask if bMVA is FPT when parame¬ 
terized by the first lower bound B. However, we can see in the previous reduction that we obtain a graph 
with S = 2, and thus the problem is even not in XP unless P = NP. 

5 Conclusion 

In this article, we presented some negative and positive results for a multidimensional binary vector as¬ 
signment problem in the framework of parameterized complexity. Notice that neither lower bounds of 
Theorem^ nor Theorem ^ are able to rule out an algorithm running in 0*(2^) (when n is part of the 
input), hence the existence of such an algorithm seems a challenging open problem. Another interesting 
question concerns the improvement of the 0{k^m) kernel of Theorem by getting rid of the parameter 
m: does bMVA admit a polynomial kernel when parameterized by k only? 
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